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ABSTRACT 

We investigate how long wavelength inflationary fluctuations can cause the background field 
to deviate from classical dynamics. For generic potentials, we show that, in the Hartree 
approximation, the long wavelength dynamics can be encapsulated by a two-field model op¬ 
erating in an effective potential. The latter is given by a simple Gaussian integral transforma¬ 
tion of the original inflationary potential. We use this new expression to study backreaction 
effects in quadratic, hilltop, flattened, and axion monodromy potentials. We find that the 
net result of the altered dynamics is to slightly modify the spectral tilt, drastically decrease 
the tensor-to-scalar ratio, and to effectively smooth over any features of the potential, with 
the size of these deviations set by the initial value of power in large scale modes and the 
shape of the potential during the entire evolution. 
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1 Introduction 


Very recent observations PEI 13 S] are on the threshold of making an important distinction 
about the details of the inflationary regime that gave birth to our universe, by the precise 
measurement of the parameter n s — 1, the slight scale dependence of the primordial pertur¬ 
bations, and the constraint on r, the ratio of the power of tensor-to-scalar fluctuations. The 
distinction arises because of the slow-roll relations 

V 72 V" V' 2 

n s = - 3 ^f + 2 — ~ °-97, r = 8—< 0.11 (1) 

written in term of derivatives of the potential driving inflation, and with observed values 
given by PEIP- From this we see that the threshold value r < 8/3(1 — n s ) ~ 0.08 implies 
that the inflaton field has a negative mass squared term, and is thus tachyonic. It is still 
too early to definitively tell whether this is the case or not, but experiments planned for 
the immediate future PM will be able to distinguish this. If the inflaton does indeed 
turn out to be tachyonic, it is worth asking what further consequences this could have for 
the dynamics of the early universe. The hallmark feature of tachyonic fields is that long 
wavelength modes are unstable (over and above the typical instabilities of light fields in an 
expanding universe). The possibility that these long wavelengths can have significant enough 
buildup to heavily influence the dynamics of the background mode was explored in [9] jlOj, 
where the parallel with a similar phenomenon in condensed matter, the spinodal regime , was 
drawn. This is the tachyonic regime, where long wavelength modes experience an instability 
that leads to exponential growth. In these and in a subsequent work m it was found 
that backreaction can be significant enough to even drive substantial periods of inflation 
on steeply tachyonic potentials, where the slow-roll conditions do not hold, such as natural 
inflation mmm with subplanckian values of the axion decay constant. Admittedly, the 
region of initial conditions suitable for spinodal effects seems quite tuned. However, if the 
consequences can be so drastic in this special case, it is worth considering whether smaller 
effects may be present in other models, especially given that data may soon definitively show 
that the inflaton is tachyonic. 

To this end, we explore whether a spinodal buildup can affect the dynamics of some of the 
typical models of inflation, and find several general features. The magnitude of the backre¬ 
action effect generically depends on the initial amount of power in the quantum fluctuations. 
Accounting for these fluctuations in the energy density does not, in general, cause them to 
exponentially dilute away, as in the standard scenario where their effect is neglected. This is 
actually a consequence of the slow-roll shape of the inflaton potential, which translates into 
the effective potential governing the strength of the fluctuations being almost equally as flat. 
If the inflaton field is tachyonic, exponential buildup of long wavelength modes does occur, 
but the effects can be seen even in nontachyonic models of inflation if the initial strength is 
great enough. We observe several effects: 

• There can be stages where the long wavelength fluctuations dominate the contribution 
to the power spectrum of density perturbations. If this is the case, the amplitude is 
not set by the overall mass scale of the potential, but instead by the initial conditions 
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for the long wavelength modes. Additionally, this shifts the value of n s — 1 by an 
amount that depends on the potential. Typically, this is a few percent for large held 
models, while in small held models the tilt is driven to be exactly scale invariant. As 
we will see, this will be enough to alter the predictions of some models to lie within 
the preferred contour in the n s — r plane. 

• Though the value of the tilt is only slightly altered, the prediction for the tensor-to- 
scalar ratio can drastically change if these effects are at play. This is because the 
scalar power spectrum can be much larger than the regular contribution, while the 
tensor power remains unchanged. This allows for inflation to proceed at a much lower 
energy scale, which decreases the tensor fluctuations. 

• The spinodal regime will be an inevitable attractor for purely tachyonic models, though 
approach may be very slow. In these models, it is possible to differentiate whether 
inflation lasted only several e-folds longer than what we observe, or if there were a 
parametrically longer phase preceding this, by noting whether the spinodal regime has 
been reached. 

• A final effect is that, if features are present in the potential, the contribution of the 
fluctuations can serve to wash them out, providing a sort of “low-pass filter” in which 
the inflaton is influenced only by the large scale trend of the potential. This can even 
make slow-roll inflation possible on otherwise completely unsuitable potentials with 
many false vacua and steep peaks. 

The spinodal phenomenon is in actuality nonperturbative, and so to make any headway 
the Hartree approximation [TS] is employed, which is a self-consistent loop resummation 
scheme. The end result of the Hartree approximation is to replace the original model of 
inflation with a two-held model, where the second held represents the power in the long 
wavelength modes. We refer to this second held as the diakyon , after the Greek word for 
fluctuation, and study its influence on the behavior of inflation as a function of its initial 
value. Even though the background dynamics is effectively two-held inflation, isocurvature 
perturbations are not generated, because at the basic level the second held decays exactly 
the same way as the inhaton during (p)reheating, so the dominant effects of this other held 
is to alter the trajectory of the inhaton, and to set the value of density perturbations. What 
we find is that in spite of the presence of drastically different regimes, where the diakyon 
held can provide the dominant contribution to the energy density, the perhaps most salient 
feature of inflation, the production of a nearly scale invariant spectrum of perturbations, 
remains robust. We expect the alterations we do observe, a drastic decrease in r and a 
slight shift in n s , to remain true beyond the Hartree approximation, at the very least at the 
qualitative level. 

Our paper is organized as follows: we devote section [2] to the technical framework. We 
begin by using the path integral formalism to derive the leading quantum corrections to the 
equations of motion. Afterwards, we discuss the Hartree approximation, how the dynamics 
is equivalent to replacing long wavelength modes with a second classical held, and show how 
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both can be easily implemented by a simple integral transform of the original potential of 
the theory. The section ends with a discussion of density perturbations in this framework. 
In section [3] we apply these results to some specific models of inflation. We begin with the 
simplest model of inflation, m 2 (f) 2 , and show that for certain initial values of the diakyon, 
the tensor-to-scalar ratio can be made much smaller without altering the other predictions 
of the model. We go on in section T2 to study the simplest model of a tachyonic held, 
hilltop inflation, and show that there are several regimes alternative to the standard slow- 
roll behavior, depending sensitively on the initial value of the diakyon. In section |3.3 


we 


look at battened potentials (0 2 / 3 and, briefly, Starobinsky), and show that the tachyonic 
nature of the inhaton necessarily causes deviations from standard inflation if the inhaton 
started from the eternal regime. I 11 section 3.4 we discuss spinodal effects in monodromy 
inflation. We find that for small values of the axion decay constant /, these serve to taper 
any oscillations in the spectrum. We also use this to highlight how spinodal effects can allow 
for slow-roll inflation even when the potential has many pocket false minima, and discuss 
this general phenomenon. We conclude in section [4} 


2 Hartree Approximation 

We devote this section to explaining the Hartree approximation as a tool in held theory, and 
justifying its use in the case at hand. We provide a simple path integral derivation of the 
Hartree approximation, which allows us to formally assess the validity of the approximation 
and, if desired, compute systematic corrections to it. Though undertaking these tasks is 
beyond the scope of the present paper, we plan to return to these in future publications. 
Afterwards, we distill the Hartree approximation down to a simple integral transform of the 
potential, allowing for efficient implementation of spinodal effects. 

A useful tool in the study of held theories is the background held method. I 11 this 
method all helds are divided into background quantities, that may or may not depend on 
space and time, and quantum fluctuations on top of these. One is then typically interested 
in two things: the evolution of the background held, and correlations between perturbations 
at different spacetime points. This method is routinely applied in inhationary cosmology, 
where an inhaton held acquires a time-dependent expectation value driving the expansion of 
the universe, and the fluctuations source perturbations in the energy density that eventually 
collapse to form galaxies. To this end, we assume the dynamics of the early universe is 
governed by a single held rolling down a potential, and perform the background splitting 

<j>(t,x) = <j>{t)+il>{t,x). (2) 

This splitting is made unique by the condition that the tadpole of the perturbation vanish, 

M = 0 =► (0) = l (3) 

To lowest order in 'll), this condition implies the classical equations of motion for the back- 
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ground field, and a linear equation determining the fluctuation: 

9, 2 0 + 3 Hd,4, + V'($> = 0, 

(_□ + V »$)) = 0. (4) 

The fluctuations are usually very small, and so then this is a good approximation. During 
inflation, however, even though the fluctuations are small, once a mode of a given wavelength 
gets stretched to larger than the horizon size it freezes to a constant value. If this occurs 
for very long time, a large number of modes will exit the horizon and contribute an effective 
stochastic offset to the equations of motion for the background, that can be encapsulated in 
the addition of a noise term DUE], the backreaction of which was recently considered in [18]. 
This can lead to drastic inhomogeneities on the largest of scales, and produces divergences 
if inflation lasts for an exponentially long time (for a review see [IS]), but for the purposes 
of a single observer measuring local correlators, this only contributes as a shift in the local 
time coordinate, as set by the value of the background part of the held. In contrast, we are 
considering the influence of long wavelength modes on the dynamics of these equations that 
does not manifest itself as a source term. 

There is a further reason to expect deviations from the lowest order during inflation if 
the held is tachyonic. In this case wavelengths larger than the inverse mass scale of the held 
will grow exponentially, which can compensate or even overcome the exponential dilution, 
and can cause significant backreaction on the value of the background held. In this case 
the tadpole condition does not reduce to the classical equations of motion, but instead 
becomes a system depending on an infinite number of the correlators of the huctuations. 
Likewise, the equations determining the values for each of these correlators will depend on 
a number of the others, as well as the background. There is no hope of solving this infinite 
set of equations for an infinite amount of variables analytically. Fortunately, there is a well- 
motivated approximation, devised by Hartree, which allows us to truncate this system to 
a manageable level, and typically yields results that are quite close to the full calculation 
(see for instance [13]). The approximation consists of treating all propagators as if they 
were moving through an external medium set by the effect of the propagators themselves. 
In this way, we acknowledge that the fluctuations of the fields can yield a large effect on the 
quantities we are trying to compute. As shown in Fig [TJ this method entails replacing all 
propagators by the sum of insertions of virtual loops, with couplings dictated by the bare 
Lagrangian. 

Once this replacement has been made, self-consistency demands that we do the same 
procedure for propagators flowing through the loops, and so on and so forth, ad infinitum. 
It should not be difficult to convince the reader that this quickly becomes impossibly hard 
to manage. Fortunately, there is a more operationally friendly method of enforcing this 
approximation, as espoused in one need only replace all higher n-point functions with 
products of propagators, multiplying two or fewer powers of the fields. For example, -0 2n —y 
Co('0 2 ) n + Ci( , 0 2 ) n-1 '0 2 j where the coefficients are combinatoric factors. Similarly, for an 
odd number of the fields if } 2n+1 —> C2('0 2 ) n, 0. The coefficients will not be important for our 
purposes, but can be found in [TO] - In this way, we have eliminated the need for correlators 


4 


0 

i i i 

1 + i + i — o 

i i i 


0 + W + oOo + " ■ 

Figure 1: Graphical depiction of the Hartree resummation. Here the dashed lines correspond 
to the background field, and the solid lines to the fluctuations. The top row is the tadpole 
condition, beginning with a line ending on a x representing the classical equations of motion. 
In the last row, self-consistency demands that once the propagator is substituted with this 
series, one must iterate the procedure, replacing every propagator in the series itself with the 
full series of propagators, and so on. The resulting diagrams lead to this sometimes being 
called the cactus approximation. 

of much higher powers of the fields, and reduced the system to a set of two equations 
determining (j) and ('ll* 2 ): 



+ 3 ffiw+f; i (Tfk ) k = o, 


k =0 


2 

2 \ \ k 


D+E W<^>) 

k =0 ' ' 


-0 = 0 . 


(5) 


These expressions were derived in [9] and their dynamics for natural inflation explored in 
[10]. A few brief comments before we go on to derive this approximation from a path integral 
point of view: the Hartree approximation is concerned with allowing the form of propagators 
to be altered to account for the background shift. In principle, higher order correlators will 
also be altered, for example the value of the quartic coupling may receive corrections by 
inserting interaction loops through different vertices. This effect is deemed unimportant 
for the study of two-point functions. Secondly, it can be seen in Fig. [l]that only contact 
insertions are included in the sum, in the sense that the (dressed) internal line originates and 
terminates at the same place on the propagator, with no further interactions on its journey. 
This ignores substructure wherein a second loop might split off from a first and recombine on 
an entirely different point. Third, if the Lagrangian has interactions of infinitely high order, 
such as the ones we will consider throughout this paper, then the Hartree approximation 
includes arbitrarily high loop order. This makes it a bona fide resummation scheme, and not 
a usual loop expansion. The path integral derivation should dispel any misgiving of giving 
precedence to extremely high order contact interactions over relatively low order noncontact 
interactions by showing how this arises naturally from a single insertion of the interaction 
Hamiltonian. 
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2.1 Path Integral Derivation 

Now we turn to a derivation of the Hartree approximation in terms of the path integral. 
This will provide the crucial insight we rely on in implementing the Hartree approximation 
in example potentials. We start by computing the tadpole of the fluctuation in the path 
integral formalism 

(ip x ) = f D ^e iS 0+^ x (6) 

The action can then be expanded in powers of the fluctuation S[(p + ip] = S + S^ip + 
S^ip 2 + Si n t, where, since we take the kinetic terms to be quadratic in fields, Si nt = 
— j d 4 x Y ^™=3 Here superscripts represent derivatives with respect to the fields, 

and barred quantities are evaluated on the background. Then 

{ip x ) = e iS J Dipe iSWi * iSi * t ip x e iSW 'P (7) 

So far this expression is exact: we have just factored out the quadratic part of the action 
in anticipation of approximating this as a Gaussian integral. Now we expand the first 
exponential 

(ip x ) « ie iS J DiP{S^'iP + S mt )' l p x e i ~ s{2) ^ 

= ^ JJ A - E j2kTTy 9t2t+1> < +1 ) (8) 

We have used that the integral of an odd number of ip' s vanishes. If we were working to 
lowest order we would also neglect the interaction term, and this would be equivalent to 
the classical equations of motion. Now we see the essence of the Hartree approximation: 
we insert a single interaction potential in the correlators we wish to compute. Higher order 
corrections to this will simply be the neglected terms in the expansion of the exponentials. 
We then treat the integrals over the Gaussian as occurring in the free field vacuum, 

J Di>'ip x 'ip y e lSW ' p2 = (ip x ip y ) o = G xy , (ip x ipy k+1 ) o = ^-^j^GxyGyy (9) 

The coincident limit of the Green’s function will be divergent as always, so we must regularize 
to arrive at finite answers. We make use of results m showing that the renormalization 
procedure in this setting does indeed result in replacing bare couplings with renormalized 
values. Our investigation will in any case be of the infrared effects of these fields, and 
so we will not be concerned with the ill-behaved ultraviolet properties of these quantities. 
Notice that we have used Wick’s theorem to express n-point functions in the free vacuum in 
terms of two-point functions. This is the crucial fact that makes the Hartree approximation 
equivalent to summing cactus diagrams. Enforcing the tadpole condition finally leads to 
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This implies that the quantity in parentheses vanishes, equivalent to the first equation in 
(§- So far we have only shown half of the equations in the path integral formulation. Before 
we discuss the other half we discuss how to replace the coincident limit of the correlators 
that appear in our expressions with a classical held. 

2.2 The Diakyon Field 

From ([5]) it can be seen that the background dynamics of the inhaton can be interpreted as 
if the held were evolving in a two-held model in the transformed potential 

00 1 / 2 \ fc 

V (> j ■ in) 

k =0 k / 

Here the second held a(y) = \jG yy represents the power in the fluctuation held. This 
second held takes on the role of encapsulating the average amount of power stored in the 
fluctuations, and has dynamics in its own right. Hereafter, we refer to it as the diakyon after 
the Greek word for fluctuation, dLanviiavarj. Thus, a appears as a second classical held in 
our equations of motion, rendering our potential effectively dependent on two helds. 

We hrst must establish the legitimacy of grouping long wavelength modes into a classical 
held such as the diakyon, as a suitable measure of the effect of tachyonic modes during 
inflation. Such a two-held interpretation was initially proposed in [j9], in which the authors 
separated the contribution of the correlator into two, based on whether the mode q is inside 
or outside the horizon. Then 

& = J(^m q < aH , ( 12 ) 

corresponding to the diakyon, was defined in such a way that it only comprises of modes for 
which q < aH. Hence, calling upon the equation of motion of the diakyon means setting 
an increasing dynamical cutoff in g-space, q max -i as a function of time, the horizon ( aH ) _1 
decreasing during inflation. A numerical analysis involving a hxed q m ax and mode equations 
(see for instance HU), on the other hand, means considering the presence of both subhorizon 
and superhorizon modes, regardless of their tachyonic characteristic. 

In HU, the choice was made not to call upon a condensation of all the unstable modes in 
a classical diakyon-like field, in favor of a cutoff in g-space. Using mode equations, however, 
one needs to make choices in order to comply with numerical limitations while still accounting 
for most of the effect due to the presence of tachyonic modes. One choice corresponds to 
qmax- Ideally, one should let q run to infinity. The validity of effective field theory (EFT) 
allows us to have a finite cutoff in q. However, the higher the q max the greater the total 
number of equations in the system, the latter itself depending also on the step size in q of 
our numerical analysis. As time goes on, more and more modes will exit the horizon and have 
some participation in the overall spinodal effect. Hence, most of the impact attributed to the 
spinodal instabilities happens during early times, as a result of the long wavelength modes 
crossing the horizon P|. At very early times, however, aH is vanishingly low, which in turn 
implies that the lowest few g-modes are responsible for the major part of the spinodal effect. 
Therefore, for such times, requiring a high q max , while seemingly more precise, may also be 
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unnecessary. Encompassing the lowest g-modes into a classical field thus looks efficient to 
incorporate most of the consequences resulting from tachyonic instabilities. 

Our goal in this paper is not to give a very precise quantitative picture but a more qual¬ 
itative one. After all, we are always limited by the accuracy of numerical computations. In 
[21] a lower bound on the total number (beyond 360 depending on the mass of the field) of 
e-folds necessary to study spinodal effects after grouping the mean field and quantum fluc¬ 
tuations into one classical field, was given, in the large N approximation. An argument was 
also made as to why such a classicalization is valid after only a few e-folds. The observation 
that the long wavelength modes are responsible for most of the spinodal effect [9], implies 
that grouping the modes which exit the horizon the earliest into one classical mode may be a 
valid approximation, especially if the total number of e-folds in the model is large. Allowing 
initial conditions that provide enough inflationary e-folds, and considering times 50 to 60 
e-folds before the end of inflation should be sufficient to give a consistent qualitative picture 
of the effects of spinodal instabilities on inflation. 

One should also address the choice of initial conditions. In refs. [9] and [21] an effective 
initial condition for the diakyon was given, 


H 

a ° = 27T 

corresponding to a zero-temperature vacuum initial state. A vacuum state is a simple but 
nonetheless valid choice of initial state. Given the previous discussion, how much the spinodal 
instabilities manifest themselves depends on the initial value for the diakyon field compared 
to the inflaton. Thus, the number of initial conditions is the same whether we choose to 
call upon the mode equations, or a two-field picture. In the former framework, one needs 
to choose q max and (p 0 while in the latter cr 0 and 0 O need to be set. Another parameter that 
needs to be fixed in both scenarios is the scale of inflation. If the scale of inflation is low, 
this allows the amplitude of (-0 2 ) to be large enough that it has a significant contribution 
to the overall dynamics of our system, while still maintaining the correct predictions for the 
other inflationary observables. 

Having discussed the regime of validity of our classical two-field framework, we can now 
proceed to discussing how to transform our potential from a function only involving one 
field, to one suitable within a two-field construction. 

2.3 The Hartree Transform 

Eq. ( [TT| is a cumbersome formula that involves Taylor expanding the potential of the field, 
then replacing each factorial coefficient by a different number. This is not only laborsome, 
but can quickly become impossible, even for some of the relatively simple potentials we 
consider here. Additionally, the convergence properties are worsened by this substitution; 
even for the simple example of noninteger monomials this transformed expression fails to 
converge. 

Fortunately, the expression can be reproduced by a very simple integral transform 
of the potential, circumventing the need to ever Taylor expand in the first place. We define 



the Hartree Transforn 0 of a function to be 


U[V]((j),a) = V u {(j),a) = / + 


y/^TTi 


na 


(14) 


It can be easily seen that the result of the Gaussian integration enforces the correct alteration 
of the coefficients to reproduce the Taylor-expanded result 0- This transformation also 
has a very natural interpretation from the path integral perspective, as it exactly coincides 
with the approximation we made there to insert a single power of the interaction potential 
in the free-held vacuum (the quadratic term, though not a part of the interaction potential, 
automatically assembles with the rest of the terms to yield the unsplit potential). In fact, we 
have the general formula: = (/('?/’))o, and can make use of the fact that derivatives 

with respect to background quantities can be pulled out of the brackets to directly rewrite all 
results in terms of the Hartree transform. The integration calls for the function to be defined 
everywhere on the real axis, which may not be part of the dynamically accessible regime, as 
in fractional powers of fields or logarithms. In order to ensure the reality of the result, we 
choose the nonunique but physically reasonable prescription that the function vanishes in 
inaccessible domains. 

Let us interpret this expression to get an intuitive understanding of what the Hartree 
approximation actually does. Here we see that we are smearing over the fluctuations with a 
Gaussian filter, with a variance that is allowed to depend on time as determined by the new 
equations of motion. This will impose that the value of the new held will be given exactly 
by the variance of the fluctuations, and will treat the evolution of the background held as 
occurring in the “homogenized” replacement held, the diakyon. Statistically speaking, this 
will give the same results as if we were to calculate the background and mode functions 
of the fluctuations, and compute the two-point corrlcator for those (to the extent that the 
approximation holds). Additionally, in the limit that the fluctuations do not have a strong 
inhuence on the background mode, the held a approaches zero, and 


H[H](0,O) = [ d#(V0 V$ + ip) = V(fi). 

Jr 


(15) 


The distributional nature of the small cr limit will make this formula cumbersome for some 
purposes, but the expansion .0 is quite easy to handle. For instance, we see that the mass 
of the diakyon at 0 is m 2 = V^, which encapsulates the fact that the strength of fluctuations 
will only grow if the inhaton is tachyonic. In this sense, applying the Hartree approximation 
will automatically tell us if the quantum corrections to the equations of motion are important 
or not. If we find that the diakyon is driven to 0 (or even the small value H/(2tt), as the 
case would be in de Sitter space) for any initial value on a given potential, we can conclude 
that the usual lowest order approximation is sufficient, while if we uncover that the diakyon 
is not pushed to small values we conclude that these effects are important. For the slow-roll 
potentials we consider, the mass of the inflaton is much smaller than the Hubble rate, which 
implies that the potential in the fluctuation direction is similarly flat. This means that even 

This is actually known as a Weierstrass transform, itself just a convolution against a Gaussian 
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for positive mass squared the field is not driven to its standard value a ~ H/(2ir) very 
quickly; in fact, for chaotic initial conditions we would expect the two values of the fields 
to be comparable, meaning that the fluctuations would not necessarily reach the attractor 
value until near the end of inflation. This is in stark contrast to the scenario in which the 
backreaction is not taken into account, where any initial power is exponentially diluted away 
and the field equilibrates to the de Sitter temperature. This translates the question of the 
importance of quantum corrections to the dynamics to a problem of initial conditions for the 
perturbations. If we begin inflation on a potential that is otherwise suitable for slow-roll, 
but the initial value of the diakyon is very large, it will continue to play an important role 
throughout the evolution. In addition, if we envision that inflation started with a fast-roll 
phase [22], in which the inflaton has a large tachyonic mass, then the diakyon quickly gets 
driven away from its standard value, and so we would expect to find that the slow-roll 
evolution proceeds with a much stronger amount of fluctuations, regardless of when in the 
inflationary history the fast-roll event took place. 

We end this discussion by noting the important identity 


aSiV^a) = d 0 V^, 


CT 


(16) 


which holds at every point on the potential. This can be verified directly, but can also be 
arrived at by the observation that if we formally replace a 2 = t, </> = x then this is the heat 
equation, and (14) resembles the heat kernel. This equation will be an important tool for 
studying the behavior of the system, and we make use of it many times in the remainder 
of this paper. For example, we can use it to find an expression for the mass of the diakyon 
that is valid for any value of the fields: 


ajv® = d'iV " + a 2 div 


(17) 


This recovers that the masses of the two fields are equal for small er, but is an exact expression 
that informs us that the mass of the fluctuation can either be significantly larger or smaller 
than the inflaton for large field values. 


2.4 Diakyon Dynamics 


We have shown that the background dynamics can be encapsulated by adding an additional 
diakyon field representing the strength of the fluctuations and replacing the potential by its 
Hartree transform. We now show that these replacements perfectly capture the behavior of 
the fluctuations as well. In addition to the background dynamics, the equations governing 
the fluctuations can be derived in the path integral formalism. To do this, it is convenient 
to briefly introduce the two-particle irreducible formalism (23) 25 US], since this is specially 
designed to study the behavior of the fluctuation spectrum in the regime where quantum 
effects can become important. This simply entails adding a bilocal source for the fluctua¬ 
tions in addition to the usual source introduced when we are only interested in one-particle 
irreducible questions, as demonstrated by 


Z[jxi 


xy J 


D'l/je iS ^ + ^ +i $ dx: > x '^ x ~^ f dxdyk xy il> x ipy 


(18) 


10 




The effective action is then defined as the Legendre transform of the logarithm of the previous 
quantity with respect to both j x and k xy . Enforcing that the derivative of Z[j x ,k xy \ with 
respect to j x vanishes yields the tadpole condition (0) = 0 that we explored above. In 
addition, we impose the gap equation, namely that the derivative with respect to k xy also 
vanishes. This will yield a nonlinear equation that dictates the spectrum of fluctuations. 
In the full theory, this will again involve correlators of arbitrarily many fields, and will be 
unsolvable. We can employ the same Hartree approximation to truncate the series to a finite 
subsystem, which yields 


(a + V& a, G'£) )a xy = i6(x - y). (19) 

If the two-point function is small, this leads to the lowest order result that the two-point 
correlator is just the Green’s function of the second variation of the action. This standard 
result can be significantly altered in regimes where the two-point function is much larger, as 
explored in [26] . where it was shown that self interactions serve to taper secular growth of 
correlators of massless fields (also see |27| . where a refinement on the Hartree approximation 
was recently employed). Since the Green’s function is just the time ordered product of the 
mode functions, G xy = > this equation is equivalent to the second equation in (J5]) . 

To show that the dynamics of these fluctuations can be replaced with the classical evo¬ 
lution of the diakyon field, we simply use ( [l6| ) to rewrite ([5]) as 

Qt + d ff P H (0, a) = 0. (20) 

With this, we have shown that the entire background dynamics can be reproduced as two 
classical fields in the altered potential H h (0, a), and have completely recovered the dynam¬ 
ics laid out in pnnQ Since we have justified to call upon a to encompass the effect of 
fluctuations, we will omit the bar in 0 and refer to 0 as the inflaton, in the rest of the paper. 


2.5 Density Perturbations 

Though the background behavior is encapsulated as a simple two-field model, the diakyon 
is in actuality a composite of the inflaton field, and so we would expect fluctuations in the 
energy density to be different from the usual two-field formula. The expression was derived 
in [10] by calculating the perturbation of the time-time component of the stress-energy tensor 
before resumming the modes into the mean field a. This is the most convincing procedure 
for obtaining the spectrum, but here we content ourselves with reproducing their results by 
a simpler method, where we employ immediate use of the diakyon approximation. Then 


- j (e[H 2 ] - 2VB.[V] + H 2 ) (21) 

2 There is also the Friedmann equation determining H in terms of the background and fluctuations. This, 
and any other equation for fields not directly involving 0, are trivially reproduced by the Hartree transform. 


1 ((H(0 + 0)-H(0)) 2 )o / 

6 H 2 \(^ H ) 2 + ( v; H 
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This expression is fully general. In the limit of small fluctuations, which is the case for 
instance if they are set by the typical de Sitter fluctuations, the perturbed potential can be 
Taylor expanded and the expression becomes 


Pk 


v- 


(vff + (v?y) 


( 22 ) 


Using an identity similar to (16) this can be rewritten as 

2 

/ U H \ 

Pk 


’* “M (yvy + (ymy I (0‘ 2 + 0‘ 4 3|)H[y 2 ], 


(23) 


which, upon taking the limit a — > 0 agrees with the result derived in m- 

The a —> 0 limit in this expression appears somewhat pathological, due to the apparent 
vanishing of the power spectrum. In this limit the standard scenario is in fact recovered, as 
the power of the large-wavelength modes can never go below the de Sitter value, meaning 
that er — y H/(2 tt), where it tracks the value of the potential. Operationally, we make the 


replacement cr 2 —» a 2 + U/(127r 2 ) in all our expressions, as justified by the definition in (12) 


If we treat <r 2 /M 2 as a small parameter, all but the leading dependence can be neglected 
in the expression for the power. In this case the Hartree transform becomes the trivial 
integration against a delta function, and the expression reduces to 


n U 2 / 2 V 
P ‘ 4 ijf + l2iV 


1 


(24) 


The interpretation of this limit is clear: the V 2 /V 2 prefactor is the usual quantity that makes 
the expression gauge invariant, and the term in parenthesis is the expected value of the two- 
point function. If the initial power on large scales is negligible, then this is just given by the 
de Sitter temperature, whereas if the power on large scales is large it is given by the initial 
value. What is nontrivial about this expression is the fact that now we are working in a 
consistent framework that dictates the behavior of large scale power for any given value after 
evolution. One crucial feature to note now is that in the limit where the diakyon dominates 
the power, the power spectrum does not depend on the overall scale of the potential. This 
is in great contrast to the usual case, where the overall scale of the inflaton potential is 
commonly set to give the observed normalization. Here this condition sets the initial value 
of the diakyon field. This translates the usual practice of trying to understand the amplitude 
of perturbations in terms of the microphysics dictating the precise potential to a question 
of why the diakyon took the value that it did. This also allows a set potential to yield a 
range of values for the fluctuations, making the observed value a potentially environmental 
quantity. 

To match the value we see, we arrive at 

< 25 > 
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so that our expansion in small a/M p is indeed justified. Note, however, that even though 
the diakyon is small compared to the Planck scale, we are still free to be in the regime where 
it can be larger than the de Sitter temperature. From these considerations, we expect that 
the diakyon will have its greatest effect in small held models of inflation, where the inhaton 
is much smaller than the Planck scale as well. 

We end this section by discussing the measurable quantities related to the power spec¬ 
trum. Still in the limit a <C M p , we arrive at 


n s — 1 = —e 



+ 2ij 


1 

1 + y 


->• 


—6e + 2 rj y —>• 0 
—4e y —> oo 


(26) 


Where we have introduced y = 127r 2 cr 2 /W as the dimensionless measure of fluctuations on 
large scales. In the large diakyon limit, the contribution from rj completely drops out! This 
tells that for small held models, where e -C 1 [28], the spectrum becomes practically scale 
invariant in this limit. For large held models, whether the tilt increases or decreases depends 
on whether e or rj is larger. In the eventuality e = rj, the prediction for the tilt remains 
constant throughout the transition to the spinodal regime. This is the case for quadratic 


potentials, and we explicitly verify this in section 3.1 


Similarly, we can derive an expression for the tensor-to-scalar ratio, where we simply 
use the fact that the value of tensor fluctuations are still set by the (Hartree transformed) 
potential 


r = 


16e 

1 + y 


—>■ 


16e 

0 


y-> 0 

y —» oo 


(27) 


which displays the alluded tendency to vanish for large values of the diakyon. Taken in 
combination, these two observables indicate that a large value of the diakyon tends to push 
predictions towards the bottom right corner of the n s — r plot. Simple application of this 
formula presumes that the values of e and rj remain unchanged no matter the value of the 
diakyon, in particular that the value of the inflaton field 50 or 60 e-folds before the end of 
inflation is unaltered. This is certainly true for small enough values of the diakyon, but away 
from this limit the field trajectories can be solved numerically. 

From these expressions we notice a general tendency of spinodal effects to push infla¬ 
tionary predictions down and to the right in the n s — r plane. This motivates these effects 
to be considered especially in models that traditionally lie to the left of the Planck con¬ 
tour, in the hopes that corrections might bring them back into agreement with observation. 
We choose to illustrate this in the simplest scenario of hilltop inflation in section |3.2[ but 


some other potentials this could potentially be relevant for include inflection-point inflation, 
Coleman-Weinberg inflation, double well inflation, and pseudonatural inflation 


3 Predictions for Specific Potentials 


We now apply our findings in the previous section to some specific potentials, to see what 
effects the backreaction might have. We begin with the simplest model of chaotic inflation, 
a single field in a quadratic potential, in section |3.1[ The Hartree transform of this model 
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consists of two noninteracting fields with equal mass. The only influence of the diakyon field 
in this setup is through the fluctuations, and can serve to drastically decrease the tensor- 
to-scalar ratio, while leaving the spectral tilt fixed. We then turn to the simplest tachyonic 
model, quadratic hilltop inflation, in section |3.2| We show in this model that the tilt is 


driven to be exactly scale invariant and the tensor perturbations unobservable, but that it 
takes typically ~ 100 e-folds to reach this regime, providing a window agreeing with the 
observed values that could possibly encompass cosmic microwave background (CMB) scales. 
A tachyonic large field model is considered in section 343 where we consider a flattened po¬ 
tential. This model eventually attains a spinodal regime for practically all initial conditions, 
allowing us to infer the difference between inflation descending from the eternal regime and 
inflation that lasted only marginally longer than the observed amount. These findings are 
summarized in Fig|2j As a final application, we study monodromy models, where we see 



Figure 2: The n s — r plot for the models we consider in this paper, along with the Planck la 
and 2a contours. Allowing for strong backreaction can serve to place both m 2 </> 2 and hilltop 
inflation predictions firmly in the preferred region, and kicks the flattened model out. This 
highlights the generic trend that tensor modes become more negligible the stronger the back 
reaction is. 


that backreaction effects serve to dampen the oscillations of the potential for small enough 
values of the axion decay constant /. This can even provide a route for the inflaton to escape 
false minima, allowing for slow-roll even in potentials that would not normally allow it. 

Before considering specific models, we first exhibit a useful general result, applicable if 
the value of the diakyon remains small compared to the Planck scale. We remind the reader 
that we only consider this regime anyway, since otherwise (21) would imply that density 
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perturbations are well above unity. In this regime, the field dynamics, in the slow-roll 
regime, becomes simply 




(28) 


Here N is the number of e-folds. Several identities can be used to rewrite the diakyon 


equation in a simpler form. For instance, the heat equation identity (16) can be used to 
express this as a linear equation in a 




(29) 


and to lowest order V is the original, a independent potential. Additionally, if the inflaton 
evolution is monotonic, we can express this as a(<p(N)) to arrive at 


da 

dtp 


14 


= -rr-cr. 
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(30) 


which has the solution 


a(N) = a 0 e 




Ivf N) d ( lo s4) 

^ = a 0 e 




= <hr 


<K N ) 


v,: 


<t> 0 


(31) 


This allows us to relate the total amount of buildup of large scale power to properties of the 
potential evaluated at specific points. It is remarkable that this equation is integrable, and 
so we reflect on why indeed this is the case. In general, one expects an exponential buildup of 


modes for tachyonic masses to be given by the integral of (30), and modes to subside during 


periods of positive mass squared. If this were all the information we had about the system, 
it would not be integrable, as there is no way to relate the amount of time spent with any 
given value of the mass to the endpoints of the evolution. However, since the dynamics is 
tied to the evolution of the inflaton, periods of large mass squared are passed through more 
quickly than periods of small mass squared, and the entire evolution can be encapsulated 
in the difference between beginning and end points. We note the pathological limits of this 
expression, where the slope of the potential vanishes: if it vanishes initially, inflation lasts 
an infinitely long time, and the diakyon attains infinite values. Similarly, if we evolve this 
to the point where the slope vanishes at the end, the diakyon subsides completely. 


3.1 m 2 0 2 

We begin with a study of the simplest model of inflation, where the potential is just a 
quadratic function of a single inflaton field to show that, even though the field is not tachy¬ 
onic, the fluctuations can have a large effect on the background dynamics. We stress that 
large values of the diakyon are not evolved to dynamically, so in order for effects to be im¬ 
portant in this scenario the diakyon must start out very large. Thus, we cannot claim the 
effects we find to be generic predictions for the model as we can in the following sections, 
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but we include this case anyway because its simplicity enables a clear demonstration of how 
effects can come into play. 

This potential is simple enough that the evolution of the fields can be solved exactly in 
the slow-roll approximation. The Hartree transform yields a two-field potential 

Cic(A °) = \ m2 (0 2 + , (32) 

which is simply two noninteracting fields with the same mass. Note that for the background 
dynamics it is possible to exploit the 0(2) symmetry to arrange that only one field direction 
has an initial value. The perturbation spectrum, however, does not respect this symmetry, 
and so we prefer to work with the physical fields. 

Using the slow-roll equations 


= ~ 2M >' 

we can solve these exactly 

/ 4M 2 1V 



-2 M 2 a, 

(33) 

a(N) = a 0 J 1 - 

AM 2 N 

(34) 



Both fields reach 0 at precisely the same time, a being just a rescaled copy of (ft. If we make 
the approximation that this is the point at which inflation ends, then we can use equation 
(21) to find 


-Ncmb — 


4 M 2 ' 


Pk = Ncmb^ + O 

1V1 p 
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n, - 1 = 


N, 


r = 


CMB 


4 m 2 
3n 2 ctq 


(35) 


In contrast with the usual predictions, now the amplitude of the power spectrum is set by 
the value of the diakyon field instead of the mass scale of the potential. If we insisted on 
the standard value, a = H/(2ir), the usual results would be recovered, and the value of the 
mass would be needed to be set to reproduce the observed spectrum. If the diakyon is large, 
however, normalization of the power spectrum implies 


~ VPjN^BM p 


6 x l(T 6 M p < 15M p 


(36) 


This indicates that the initial value of the inflaton is much larger than the initial value of the 
diakyon, but the diakyon still has a much larger value than the usual case for generic m 2 . 
Evidently, the tilt is still only sensitive to the total number of e-folds, and so its standard 
predictions remain unaltered. The only parameter that depends on the mass scale here is 
the tensor-to-scalar ratio, which now decreases when m 2 takes small values. This has the 
consequence that the predictions of this model can be brought into agreement with the 
preferred values of these parameters, if the mass scale is sufficiently small. The origin of this 
result resides in the fact that, if backreaction is large, fluctuations are no longer tied to the 
expansion rate, but can now be much larger. Since the strength of the scalar fluctuations is 
dictated by the value of a, whereas the tensor fluctuations are still prescribed by the overall 
energy density, the scale of the potential can be lowered arbitrarily, while still maintaining 
the observed power spectrum. 
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3.2 Hilltop Inflation 

While the quartic hilltop model fits with the measured values of inflationary parameters well 
enough to merit inclusion on the Planck 2015 n s — r plot [I], the simplest model, quadratic 
hilltop inflation, predicts too small an n s . As such, it is the ideal case to study. We find 
already in this simple scenario that there are many different regimes, corresponding to the 
initial value of the diakyon. 

The transformed potential in this case is 

Kmtop^ v) = V Q - + a 2 ) (37) 

and, similar to the previous case of m 2 0 2 , near the top of the hill the two fields behave 
independently of each other, each experiencing exponential growth: 

(f>{N) = </> 0 e |r?o|7V , a(N) = a 0 e MN , (38) 


where rj o = —m 2 M 2 /Vo is the slow-roll parameter at the top of the hill. For simplicity, we 
consider the predictions for the case when inflation ends at the point where the potential 
becomes negative. In actuality, we expect the form of the potential to be altered well before 
this, but this simplistic setting will be enough to illustrate the effects of the diakyon field. In 
this scenario the vacuum energy can be approximated by the constant term in the potential, 
with corrections becoming important only in the last e-fold or so of the evolution. The value 
of the inflaton at CMB scales is set to be 0 cmb = \j2/rjQe~\‘ n °\ NcMB . The initial value of the 
diakyon must be less than this, or else it would spoil inflation before the required amount of 
e-folds has elapsed. 

The power spectrum is given by 

p Vo 2 * V + F) 

m 4 (0 2 + a 2 ) 2 



The previous requirement that a < <f> doubly suits us, because as it can be seen from above 
it also corresponds to the regime where Pp- —y 3/(4 r/ 2 ) 1 , certainly of no phenomenological 

relevance. In the limit <r<C^ the power spectrum becomes 


Pk —*■ 


* 

m A <fr 2 


(40) 


where it is understood that for small values of the diakyon this expression should recover 
the standard result by saturating to the de Sitter temperature. This form displays the scale 
invariant behavior uncovered in the general case, owing to the fact that both fields have the 
same exponential dependence on e-fold time (making this regime also of no phenomenolog¬ 
ical relevance). If, instead, the diakyon starts at values much smaller than the de Sitter 
temperature, we recover the standard predictions 


n s — 1 = 6 ? 7 0 e 2%(iVcMB " 7V) + 27 /q, r = -16 r /o e 2 ' 7o(iVcMB " 7V) . (41) 
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The approximations we have used break down for tjq <C but in this case the amount of 

e-folds is so great that the last 60 will occur at the base of the potential, which is practically 
linear, and does not really resemble hilltop inflation anyway. The interesting regime is when 
the diakyon starts smaller than the de Sitter temperature and then, through the course of 
its exponential growth, transitions to the regime of perfect scale invariance and negligible 
tensor-to-scalar ratio. This transition actually takes place on the order of 50 to 100 e-folds, as 
shown in Fig [3j so that this is not a singular event that would require tuning to occur within 
the observable window. This also ensures that the running of the power spectrum remains 
unobservably small throughout the evolution, in agreement with current observations. 



Figure 3: Evolution of the tilt with two different values for the slow-roll parameter 770 , the 
solid curve corresponding to r] 0 = —1/25 and the dashed to r] 0 = —1/30. Different initial 
conditions and overall scale of the potential will shift the transition region to the left or 
right. These curves illustrate that the transition takes a few dozen e-folds to pass through 
the observable window, a naive theorist’s reconstruction of which, based off [IJ, is shown as 
the shaded region. 

We can comment on other versions of hilltop inflation before we conclude this section. 
For the case of quartic hilltop inflation, which lies well within the best fit region in the Planck 
data, the spinodal effects produce 

Quartic hilltop cr) = V 0 -X (</> 4 + 6 a 2 (j) 2 + 3cr 4 ) . (42) 

So that the diakyon direction is actually much steeper than the inflaton direction. This 
necessitates the initial values to be even more fine tuned for this to match observations. 

3.3 Flattened Potentials 

There are certain contexts where integrating out heavy degrees of freedom serves to flatten 
the inflaton potential from an initially steeper profile for large field values [3D]. This can 
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cause the field to be driven by an effective potential that is a fractional power of the field. 
A typical example of such a potential that will be testable in the coming years is 


VW0 = 


(43) 


The standard predictions for this potential are 


<KN) 
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(44) 


with N 0 being the beginning of inflation and Nqmb the amount of enfolds since CMB scales 
left the horizon. However, on this potential the field is tachyonic, and so we should expect a 
region of initial conditions where spinodal effects can become important. To see this we apply 
the Hartree transform to this potential, yielding a result in terms of confluent hypergeometric 
functions 




— 10 / 3 


T (5/6) cr 2 / 3 Ml 
Ml — x! 7b 


2 2 / 3 tt 1 / 2 


3’2’ 2cr 2 


2X , T (4/3) (j) „ (1 3 

+ XT la , /o , ( jh X, X) 


2 l/6 7r l/2 a l/3 J 
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(45) 

This shows the true utility of (14), allowing for an analytic expression for the full potential. 
The latter is displayed in Fig. |4j. 



Figure 4: The full potential for 0, a, overlaid with the a -independent original potential. For 
small initial values of the diakyon the inflaton dynamics exactly corresponds to the usual 
case, and the diakyon experiences an instability towards larger field values. 

The previous expression is not very inviting of an analytical treatment, but it is easy 
enough to work with numerically, and in asymptotic limits. For small values of the diakyon 
the original potential is recovered, as must be the case. The dynamics of several choices of 
initial conditions are shown in Fig. [5], 

However, for Planckian values of the diakyon, the amplitude of the power spectrum is 
well over unity, which is nowhere close to the observed value. Therefore, we restrict our 
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Figure 5: Inflationary trajectories during the last 60 e-folds for different initial field values. 
The solid lines correspond to the inflaton, and the dashed lines correspond to the diakyon. 
For small initial values of the diakyon, as in the green curves (solid top and dashed bottom), 
it never conies to dominate, and inflation proceeds as usual. The purple (intermediate) and 
orange (coincident at start) curves exhibit extreme values of the diakyon, where it comes to 
dominate at late times, leading to a few extra e-folds of inflation. These values lead to a very 
pathological amount of power, so should be treated as exaggerated scenarios demonstrating 
the effect present in smaller quantities for more modest values. 


analysis to values a -C M p < 0, for which case the potential can be approximated as 


(46) 

In the approximations we have made, the inflaton trajectory is barely altered from the 
standard result (44), and the diakyon trajectory is 

0b 


a(N) = 


1 - 

1 No) 


1 / 6 - 


(47) 


This blows up at the end of inflation, but the approximations we have made break down 
well before that, as can be seen in Fig. [5j In the regime where a 0, we can use the fact 
that i.Fi(a; 6; 0) = 0 for all a, b to find the limiting form of the potential 

rC(°.b^b^|f‘ ll,/ V /3 . (48) 

which is of the exact same form as the original potential, but with an additional numerical 
factor fa 0.4. Before this regime is reached, though, the diakyon grows with time, character¬ 
istic of tachyonic behavior, and directly opposed to the lowest order approximation, where 
it tracks the value of the potential which decreases with time. This causes an interesting 
behavior: there is a discernible difference between inflation on this potential that comes from 
an eternal regime, and inflation that started only several e-folds outside of the observable 
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window. To highlight the difference between these two scenarios, let us assume that inflation 
began at the eternal regime 1 [31], with the initial value of the diakyon set by the de 


Sitter temperature (13): 


/4\ 3/8 Mp /4 /2" 3/2 5/2 2 1 [2 /i 5 / 2 

( 9 ) /i 5 / 4 ’ V °-\/3 M r * ’ a ° “ 12^2 V 3]^p2‘ 


(49) 


For simplicity we use the approximation that the fields start on their slow-roll attractor 
behavior, with negligible initial inhomogeneities. The total amount of e-folds from this 
regime is N e = (M p / /i) 5//2 /\/6, much larger than the observed 60 or so e-folds. By the time 
the fields have evolved down to CMB scales, the value of the diakyon is 
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(50) 


The diakyon is still subplanckian, so our expansion is valid, but it is much larger than its 
vacuum value. In this regime the fluctuations behave purely spinodally. Normalization of 
the power spectrum at N CMB = 60 then requires |i= 1.6 x 10~ 5 M p , and the predictions for 
the model become 


= 1 — 


3iV C MB 
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(51) 


CMB 


which are clearly different from (44), with half the predicted tilt and practically vanishing 
tensor-to-scalar ratio. These values of the parameters actually lie substantially outside of 
the observed values, making this regime in disagreement with experiment. 

There is a general trend that is indicative of all tachyonic models with an eternal regime: if 
inflation proceeds for a sufficiently long time, the diakyon takes over. Contrary to the usual 
predictions for large field models in this class, we would not expect to see tensor modes, 
and if application of the usual expression for the tilt were applied, it would lead us to the 
wrong conclusion about the parameters and shape of the potential. If we were confident 
that inflation took place on a given potential, however, observations would allow us to infer 
whether inflation lasted parametrically longer than the observed amount by measuring the 
amount of deviations from the standard predictions. 

We would like to conclude our analysis of flattened potentials with a brief remark on 
another scenario: the Starobinsky ‘1R + .R 2 ” model (32]- This analysis is particularly relevant 
here since it is a small field model, so if spinodal effects are present the spectrum looks 
scale invariant. If we assume that the diakyon starts at the de Sitter temperature, then 


straightforward application of (31) yields 


Wot V (H \ 2 


(52) 


Hence, if inflation lasted any longer than the minimum required to solve cosmological prob¬ 
lems, the backreaction due to the diakyon would be important and engender a shift in n s 
towards exact scale invariance. 
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3.4 Monodromy Inflation 

As a final application of the Hartree transform we consider axion monodromy inflation [33l 
M\- In the simplest realization of this scenario, inflation is driven by an axion field in which 
stringy effects break the exact periodicity of the potential. This leads to an approximately 
linear potential with small modulations, 


trnono (0) = AV + A 4 COS (j'j • ( 53 ) 

Inflation is driven by the linear term, yielding standard predictions for inflationary parame¬ 
ters: n s — 1 = — 3/(2A^ CM b) and r = 4/iVcMB- In addition, the sinusoidal correction leads to 
characteristic oscillations in the power spectrum that are logarithmic in scale [35] and have 
yet to be found at a significant level [36], but can be arbitrarily small depending on the pa¬ 
rameter A. One is free to take natural values of the axion decay constant / 
in which case the field undergoes many oscillations in a single e-fold, or the cycles could be 
so large to encompass the entire range of observed CMB scales, possibly providing a running 
spectrum that explains lack of power both on large scales and in the matter power spectrum 

ra- 

Monodromy inflation is a large field model, with 

(f>(N) = \J 2(iV cmb - N)M P . (54) 


If spinodal effects are not too large, the normalization of the spectrum enforces fi — 6 x 
10~ 4 Mp. The other parameters can be expressed in a convenient combination b = A 4 /(/r 3 /). 
If this parameter is too large, b > 1, the potential will not be monotonically decreasing, but 
will instead develop pockets that, if too deep, could potentially cause the inflaton to become 
trapped in a false minimum instead of undergoing slow roll. 

We would like to investigate to what extent spinodal effects might change this story. 


Taking the Hartree transform of the potential (53) yields 


<r) = /p 3 <P + A 4 cos( 0 //)e a * /2/ \ 


(55) 


The linear term remains untouched, and the cosine, as an eigenfunction of the Hartree 
transform, is now multiplied by a gaussian. The potential is plotted in Fig ( 6 ] . 

We see that allowing for tachyonic modes to build up and backreact on the background 
evolution can serve to taper the size of the oscillations, depending on the initial conditions 
taken for the inhomogeneity. An important point is that spinodal effects prevent the usual 
restriction to an exactly homogeneous slow-roll evolution from being an attractor, so that 
only for a narrow range of field space a < / will these oscillations occur at all. For generic 
values of the diakyon, the potential will look exactly linear, with no additional signals. As 
we have argued in section [ 2 ] if inflation is to have come from an eternal regime, then we 
would expect a much larger initial value of the diakyon field. We note that spinodal effects 
can provide a route for slow-roll inflation to occur even for large values of b , which would 
have been traditionally excluded. 
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Figure 6: The two-field potential. For large values of the diakyon the oscillations are washed 
out, and the inflaton direction is strictly linear. 


The line a = 0 corresponds to the traditional slow-roll path, which gives the strongest 
possible oscillation signal. From the form of the potential one can expect that during the 
periods when the inflaton </> is near the top of a hill the fluctuations roll towards the plateau 
region, corresponding to the buildup of tachyonic-wavelcngth modes. Before this effect 
can get too strong, however, the inflaton goes through a valley region, at which point the 
diakyon tends towards decreasing values. If the starting point is well out on the tail of the 
Gaussian, the diakyon will barely change its value at all, and the evolution of the inflaton 
will look effectively linear. These qualitative effects can be encapsulated in an analytical 
approximation below. 

To get a stronger understanding of how the diakyon behaves during inflation, we make 
a few simplifying assumptions to make its time dependence analytically solvable in terms 
of special functions. We first take the time dependence of the inflaton to be given by its 
unperturbed slow-roll trajectory (54), which is valid to lowest order in slow roll and b. 
Secondly, we take a to be in slow roll as well, which is justified because the slope of the 
Gaussian is ^ < 6^, and so is comparatively shallow to the inflaton direction unless 
b 1. Third, if the Hubble rate is dominated by the linear term in the potential, then we 
can write 


|;v^o£ = £«» 
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(56) 


which we can integrate to get the following expression for a (IV) in terms of the exponential 
integral Ei(x): 


Ei 




—26sin f^^No-N)) . 


(57) 


From this expression we can remark that the field oscillates between two extremes cor¬ 
responding to the maxima and minima of the sine function, so that one full cycle brings the 
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Figure 7: Implicit solution to equation (57) for example values b = 0.5 and cr 0 = 0.7/. 
The diakyon oscillates between the two dashed curves during each cycle. Increasing the 
initial value of the diakyon field corresponds to shifting the horizontal lines upwards, while 
increasing b corresponds to increasing the vertical width between them. Both for small and 
large initial values, there is hardly any change in the field at all. 


(x-field back to where it originally started. Moreover, the amplitude of this oscillation cycle 
does not change, and there is no long-term drift of the central value. This curve, along with 
extremum values for the field, is plotted in Fig. Jr] for example values of the parameters. 

Notice that the oscillation can be quite lopsided. If we increase the value of b , corre¬ 
sponding to more pronounced hills and valleys, the range of a will be greater. For very small 
initial values, the maximum value of a is negligible, and spinodal effects can be ignored. This 
corresponds to placing the trajectory at the very top of the hill, where there is no strong 
slope to provide a force away. Very large values of the field are plotted on the right plot, for 
which the field does not move very much at all, regardless of the value of b. In this regime 


asymptotic limits can be made in (57) to find the total range of a during one cycle to be 


Act —>- bcr 0 e 


~<T 2 /2f 2 


(58) 


This corresponds to regions far enough out on the Gaussian tail that the potential looks 
exactly flat, and so the diakyon is roughly stationary. 

We now discuss the parameter values for which the diakyon can effectively screen the 
small scale perturbations. If the fluctuations in (24) are dominated by the potential and not 
the diakyon, we require a 2 <C V/M 2 , which translates into a < 4.5 x 1CT 6 M P . In the event 
that the diakyon dominates the fluctuations, which will occur if the scale /i is smaller than 
the standard value, correct normalization enforces that this equality be saturated. In this 
case the oscillating value of a can reintroduce oscillations in the power spectrum, but these 
will be large only if cr m /, which is an apparently fine-tuned scenario that we disregard. 
Screening becomes efficient for a > 4/, which means that the mechanism we have described 


can only be operational for small values of /, say / < 10~ 6 M P . This value may seem 
somewhat low from the perspective of high energy physics models, but corresponds to the 
upper value of the “classic window” for axion dark matter [58], Monodromy models with 
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different powers of the inflaton are sometimes considered [32]- Taking into account these 
models does not substantially raise the value of / needed for the diakyon’s screening of small 
scale perturbations to take place. However, the dynamics may be more interesting, as the 
mass of the inflation causes the diakyon to increase (or decrease) with the evolution, which 
could cause the oscillations to appear only at early (or late) times. 

4 Conclusions 

We have considered the influence of backreaction from inflationary fluctuations for a few typ¬ 
ical models of inflation. We showed that the strength of the influence is largely dependent on 
the initial conditions of the inflaton, and whether the dynamics are tachyonic enough to drive 
the diakyon to large held values. If the long wavelength modes are initially populated close to 
the de Sitter equilibrium values, backreaction is negligible and inflation proceeds as normal, 
unless a very long time elapses since the beginning of inflation. If the modes are initially 
populated to a much greater extent, then it can significantly alter the inhationary dynamics 
and predictions, by inhuencing the strength of the generated perturbations. Though the 
details of inhationary dynamics can be substantially modified, the basic predictions of the 
scenario, mainly a stage of near exponential expansion, generating a nearly scale invariant 
spectrum of perturbations, remains the case. We note four general trends. The predicted 
value of the spectral tilt deviates slightly from the original value in large held models, and 
moves strongly towards scale invariance in small held models. Secondly, since the size of 
the fluctuations is no longer set by the energy scale of the potential during inhation, infla¬ 
tion can proceed at a much lower scale and still produce the observed size of huctuations. 
This causes the tensor-to-scalar ratio to drastically decrease. Thirdly, the eventual tendency 
towards the tachyonic regime allows us to infer the duration of inhation before observable 
scales. Lastly, if the potential has any sort of features, the strong backreaction regime serves 
to smear out these features, providing a smoother potential for the inflaton. This can occur 
even for features so deep that they would have prevented the inflaton from reaching its true 
minimum at all. This was studied in the simple case of monodromy inhation, but holds in 
generic potentials, as displayed in Fig. |8j 

We remark that this washing out of features is not enough to completely solve the r] 
problem, which is the question of how the mass of the inflaton can be lighter than the Hubble 
scale, even after quantum corrections are taken into account. This stems from the fact that 
the size of the features that are washed out depends on the value of the diakyon. If we 
Fourier expand a potential, H(0) = f dnV n cos(n0/M p ), the Hartree transformed potential 
multiplies each coefficient by a gaussian e - n2 °"V M p ) which effectively cuts off the value of the 
integral at n ~ M p /a. If all coefficients are of the same order, then this yields 7] ~ M 3 /cr 3 , 
but since we require a < \JeP^M v we arrive at rj > 10 18 3> 1. The only way to circumvent 
this conclusion is if the coefficients of the Fourier series become tiny for suitable values of 
n, but this requires working on a potential where the 7] problem was solved originally. As 
it stands, this mechanism is only suitable for removing features that are much smaller than 
the typical held range. 
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Figure 8: Spinodal effects in a generic random potential. For small values of the diakyon the 
shown potential is completely unsuitable for slow-roll inflation, containing many false minima 
and steep peaks. Backreactions tame these features by smearing over the original potential 
with a Gaussian filter. This picks out the large scale trend of the potential, providing a route 
for slow-roll inflation to proceed. 


The complicated nature of the strong backreaction regime forced us to employ the Hartree 
method to encapsulate the qualitative aspects of these effects. We expect it to be a reasonable 
estimate of the dynamics based on its success in other scenarios. The overall effect of this 
approximation is to replace the original single field dynamics by another field that governs 
the strength of the fluctuations. It is important to make an estimate of the next order 
corrections to this approximation to ensure that the calculations we performed are under 
full control. The path integral approach to the Hartree approximation is a simple way to 
do this, and we plan to return to this in future publications. Nevertheless, we expect the 
qualitative features we have uncovered to remain true in the full theory. 

There are a number of other issues worth exploring. We employed the Hartree approx¬ 
imation in the uniform curvature gauge, where all perturbations were in the inflaton field. 
In general, there is no guarantee that resummation schemes are gauge invariant, and so it 
would be interesting to extend this framework to gauge invariant perturbation theory. The 
Hartree approximation is also known to not respect global symmetries [3D], and in the sim¬ 
plest settings provides spurious masses to Nambu-Goldstone bosons of broken symmetries. 
As the inflaton can be thought of the Nambu-Goldstone boson of time translation invariance 
HU. it might be worth considering if a modified version of the Hartree approximation such 
as used in m is more appropriate. We have also only considered strong backreaction effects 
on the scalar two-point correlator. It might be interesting to examine if its effects on the 
tensor power spectrum can be similar, and what effect it might have on nongaussianities. 
The latter consideration might even necessitate a generalization of the above scenario to 
a 3-particle irreducible formalism, similarly causing the strength of nongaussianities to be 
dependent on initial conditions. We leave these as questions for the future. 
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